Hadron Yield Correlation in Quark Combination Models in High-Energy AA Collisions 
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We study the hadron yield correlation in the combination models in high-energy heavy-ion collisions. We 
derive the relationship between the average yields of different hadrons produced in the combination of a system 
consisting of equal numbers of quarks and antiquarks. We present the results for the directly produced hadrons as 
well as those for the final hadrons including the strong and electromagnetic decay contributions. We also study 
the net quark influence by considering the case when the number of quarks is larger than that of antiquarks. We 
make comparison with the data wherever possible. 
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I. INTRODUCTION 



Hadron yield correlations, measured by the ratios of 
the average yields of different hadrons produced in high- 
energy reactions, are one kind of characteristic proper- 
ties of hadronization mechanisms. It is usually expected 
that these correlations are more or less independent of the 
particular model especially when the combination mech- 
anism is concerned. Such properties were therefore con- 
sidered as a good probe for the hadronization mechanism 
in different high-energy reactions already in the 1970s to 
1990s [1-8]. The study of these correlations has attracted 
much attention in heavy-ion collisions recently [9-15] be- 
cause they are considered as a probe to test whether the 
(re-)combination mechanism is at work. This is interest- 
ing because whether the combination mechanism works 
might be considered as one of the signatures for the for- 
mation of the bulk color-deconfined quark matter system 
before the hadronization takes place. Experimental re- 
sults are available from the Relativistic Heavy Ion Col- 
lider (RHIC) [16, 17], from relatively low-energy colli- 
sions such as those obtained by NA49, NA61, and CBM 
Collaborations at the Super Proton Synchrotron (SPS) [17- 
20], and more recently from the very high energy reactions 
at the Large Hadron Collider (LHC) [21, 22]. These re- 
sults seem to suggest a dramatic change for collisions from 
lower to higher energies, which is considered as one of the 
hints for phase transition [18]. 

In order to make a judgment whether the combination 
mechanism is at work by comparing the theoretical re- 
sults with the corresponding experimental data, it is im- 
portant to see whether, if so, to what extent, the theoreti- 
cal results depend on the particular model used in obtain- 
ing these results. There are many studies that have been 
made in the literature [9-13, 23]. These studies are usu- 
ally based on some particular (re-)combination or coales- 
cence models and/or some particular assumptions. It is not 
clear whether the results obtained depend on the particular 
assumption(s) made in these particular models. For this 
purpose, in this paper, we will make a systematic study of 
the average yields of different identified hadrons and their 
relationships obtained in the combination mechanism. We 
will start the study by considering the case for the combi- 



nation of a system of quarks and antiquarks from the basic 
ideas of the combination mechanism. We will make the 
study as independent of the particular models as possible 
but present the assumptions and/or inputs explicitly when- 
ever necessary. 

The rest of the paper is organized as follows. In Sec. II, 
we derive the formulas for calculating the average yields 
of hadrons and their relationships in the combination of 
a system of quarks and antiquarks. We consider a sys- 
tem where the number of quarks and that of antiquarks are 
equal and discuss the net quark influence as well. We com- 
pare the results with the available data in Sec. III. Here, 
the experimental results from the LHC [21, 22] are taken 
as an example to test the predictions for the case where 
the net quark influence is considered as negligible, while 
those from the RHIC and the SPS [18, 24-36] are used to 
test the net quark influence. A short summary is presented 
in Sec. V. 



II. HADRON YIELD RATIOS IN THE COMBINATION 
MODELS 



In this section, we begin with the general formalism of 
hadron yields in the combination models based on the ba- 
sic ideas. For this purpose, we start with a quark-antiquark 
system as general as possible. Then we simplify the re- 
sults by using some explicit assumptions, simplifications, 
and/or approximations. We present the results for the ratios 
of the yields of hadrons directly produced as well as those 
including the contributions from the resonance decays. 



A. The general formalism 



We start with the most general case and consider a sys- 
tem of N q quarks and Nq antiquarks. We denote the num- 
ber of quarks of flavor q t by Nq i so that Y*iN qi = N q an d 
similarly = Nq. These quarks and antiquarks com- 

bine with each other to form the color singlet hadrons. The 
number density of the directly produced hadrons is given 
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by 

fMj(pMj)= 2 I dpidp2f qi q 2 (puP2;N qi ,N 9i )x 

9l 92 J 

( R M „q i Z ll (PM P P\,P2\ N qi ,N 9i ), (1) 

fBj(PBj) = 2 I dpidp 2 dp 3 X 

919293 J 

fq,q 2 c b {P\, P2, P3;N q .,Nq.) X 
^,9,9293 Pl,P2,P3\ N q . , N q ,), (2) 

where /m ; and /b. are the momentum distributions for 
the produced meson Mj and baryon Z?y, respectively; 

fqiq 2 (PUP2;N g .,Nq.) and fq l q 2 q,(puP2,P3\Nq i ,Nq i ) 

are the two- and three-particle joint momentum dis- 
tributions for (qiq 2 ) and {q\q 2 q3), respectively. The 
kernel functions KM^q^iPM,, Pi, P2',N qi ,N qi ) and 
K Bjmcnq ^{p Bj , Pu P2, P3;N q „N qi ) stand for the probability 
density for q\ and q 2 with momenta p\ and p 2 to combine 
into a meson Mj of momentum p M . and that for qi, q 2 
and q 3 with momenta pi, p 2 , and /?3 to coalescence into 
a baryon Bj of momentum p B .. Here, in the arguments, 
we use N qi and Nq. to represent the dependence of these 
functions on the numbers of the quarks and antiquarks of 
different flavors, and also on the total collision energy V* 
of AA reactions. We note in particular that not only the 
joint distributions f qiq2 and f qiqm but also the probability 
densities f RM j .q i c P an d f^Bj.qiq 2qi are m general dependent 
on N qi and N q .. This is because, for finite N m and Nq., 
the probability for a given quark q\ to combine with a 
specified antiquark q 2 to form a specified meson Mj or 
two other specified quarks q 2 q 3 to form a baryon Bj is in 
general dependent on the number N C[I of existing quarks of 
different flavors and the number N q . of antiquarks. 

We note in particular the relationship between the de- 
scription presented here and those given in the literature 
in different models based on the combination mechanism 
such as the coalescence model [9, 10], the recombination 
model [11, 12], and the quark combination model [13, 15, 
37] developed by different groups. Eqs. (1) and (2) are in- 
tended to be the general formulas based on the basic ideas 
of the combination mechanism. The different models are 
examples of the general case that we considered in these 
equations. In these models, different method(s) and/or as- 
sumption^) are usually introduced to construct the precise 
form of the kernel functions < Ru j ,q^q 2 {pM P pi, P2',Nq i ,Nq i ) 
and HBjA^qAPBp Pu P2, P3\N qi ,N qi ) in order to provide 
a good description of different properties of the hadrons, 
such as momentum distributions and so on. For example, 
in the recombination model developed by Hwa and collab- 
orators [12], these kernel functions are just the recombina- 
tion functions. 

The joint distributions f qiqi and f qi q 2 q 3 are the number 
densities that satisfy 

/ fq, qz(PUp2\ N qi , N qi )dp idp2=N qiq2 , (3) 

/ fqmqiiPu P2, P3\ N qi ,N qi )dpidp 2 dp 3 = N qiq2q „ (4) 



respectively, where N qiq2 = N qi N q2 , and 

C N qi N q2 N q} for qi +q 2 + q 3 

N qim3 = \ N qi (N qi - l)N q3 fox qi =q 2 + q 3 (5) 

{ N qi (N qi - l)(N qi - 2) for qi = q 2 = q 3 

are the numbers of all the possible (qqYs and (qqqY& in 
the bulk quark-antiquark system that we consider. For the 
convenience of comparison, we rewrite them as 

f qiqi (puP2;N qi ,N qi ) = N qiq2 tf;l_( Pl ,P2;N qi ,N qi ), (6) 
f qiq2q3 {pup2,P3;N qi ,N q .) = 

N 9l9 ^f^U(pi,P2,P3',N qi ,N 9l ), (7) 

so that the distributions are normalized to unity where we 
denote by using the superscript (n), i.e., 

/ Jtl 2 (PuP2;N qi ,N qi )d Pl dp 2 = 1, (8) 
/ Jtl 2q3 (PuP2,pr,N qi ,Nq t )d Pl dp 2 dp 3 = 1. (9) 

In terms of these normalized joint distributions, we have 
= JXfc f dpidp 2 f%l(puP2;N q „N qi ) 

9192 J 

x n Mj {p Ml ,p\, P2\ N qi ,N 9i ), (10) 
fB,(PBj) = J] N919293 I dpidp 2 dp 3 X 

919293 

Jtl 2q3 (puP2,P3;N qi ,N qi )x 

R Bj (PB P Pl,P2,P3;N qi ,N qt ). (11) 

Integrating over pu- or p B . from the momentum distribu- 
tions, we obtain the average numbers of the directly pro- 
duced mesons Mj and baryons Bj as 

N Mj (N qi ,N qi ) = ^N qiqi I dp M ,dp l dp 2 x 

9192 J 

jtl^'PZ' Nq l ,N qi )'R Mj (pMj,Pl,p2;Nq i ,N qi ), (12) 

N Bj (N qi ,N q .) = ^ N qmqi I dp Bj dpidp 2 dp 3 x 

919293 

4"? 2?3 (Pi- P2, P3\ N qi ,N q .)n Bj {p Bj ,pi, p 2 , P3\ N qi ,N q .). 

(13) 

For a reaction at a given energy, the average numbers of 
quarks, (N qi ), and those for the antiquarks, (N qi ), of dif- 
ferent flavors are fixed. The numbers of quarks and an- 
tiquarks follow a certain distribution which we denote by 
P{N qn N qi , (N qt ), (N qi )). The average yields of mesons and 
baryons are given by 

(N Mi K^)= J] P (N q , , Nqj , {N q , ), {N qi ))Nmj (N q , , Nq.), 
N q .N it 

(N Bj )( \ft ) = J] P(N qt ,N qi , (N qi ), (N qi ))N B .(N qi ,N qi ). 

N q ,N q , 
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These equations are the general formalism for calculat- 
ing the average yield of a certain sort of hadrons in high- 
energy reactions based on the basic ideas of the combina- 
tion mechanism. More specific results can be obtained for 
special cases when special assumptions are made about the 
distributions and/or the kernel functions. We present such 
cases step by step in the following. 



B. Factorization of flavor and momentum dependences 

The flavor dependence of the kernel functions < RM j ,q i q 1 
and < RB 1 ,q i q 2 q 3 is responsible for flavor conservation in the 
combination process and the differences between the com- 
bination probabilities for different flavors of quarks, anti- 
quarks, and hadrons. In general, the momentum and the 
flavor dependencies of these kernel functions are coupled 
to each other. In that case, the results for the ratios of the 
average yields of different hadrons can be dependent on 
the way of coupling. In this paper, we do not consider such 
coupling effects. In contrast, in the following, we consider 
only the simplest case where the momentum and flavor 
dependencies of the kernel functions are decoupled from 
each other. In other words, we consider the case where 
they are factorized, i.e., 

92 (PMj ,PUP2', N q , , Nq. ) = 
K% qi qS N <l>' N 9^ < J>M,Pl, PI', N q ,Nq), (14) 
KB jAl q 2 q z {pBj,PuP2, P3\N qt ,N q ^ = 

K^q^N^ipB^uPl, P3lN q , N q ), (15) 

where the flavor-independent parts 'R^ipM, p\ , pi', N q , N q ) 

and 7?^ (p B , Pi , P2, P3\ N q ,N q ) denote the probability for a 

(qq) with momenta p\ and p2 in a system consisting of N q 

quarks and N q antiquarks to combine with each other to 

form a meson M with momentum pm and that for a (qqq) 

with momenta p\, p2, and pj in the system to combine 

with each other to form a baryon B with momentum ps, 

respectively. The flavor-dependent parts K~J qiqi (N q ,,N qi ) 

( f) J 
and Kg \ 2 } (N qt , N 9i ) represent the probability for the q\ 

and q2 to combine into the specified meson Mj in the case 
that they are known to combine into a meson and that for 
the q\, qi, and q$ to combine into the specified baryon B ; 
in the case that they are known to combine into a baryon, 
respectively. They are taken as satisfying the normaliza- 
tion condition 

j 

y# =i. (i7) 

/_l Bj,qiq 2 q 3 V > 

j 

We further assume that the normalized joint momentum 
distributions of the quarks and/or antiquarks are flavor in- 



dependent, i.e., 

f!H(PuP2;N q „N qi ) = 4f(puP2;N q ,N q ), (18) 
f£l*<Pi.P2,P3;N 9 ,N*)=fM(puP2,P3lN q ,N 9 y (19) 
Under these two approximations, we have 

N Mj (N q „N qt ) = ^N qiq2 K < fi ]qih (N qi ,N qi ) x 

J dp M dp i dp 2 jf q ) (pi,p2;N q ,N q )X 

K%\pM,Pi,P2;N q ,N q ), (20) 

N B .(N qi ,N qi )= £ N qiqiq3 ^ q2q3 (N qi ,N qi )X 

J dp B dpidp 2 dp3f!£ > q (pup 2 ,P3',Nq,N q ) x 
■R ( £\pB,PuP2,py,N q ,N q ). (21) 

We denote 

jM(Nq,Nq, V?) = J dp M dp\dp2 x 

4f(Pi , P2; N q , N q )n^\p M ,px, P2; Nq, N q ), (22) 
y B (N q ,N q , V*) = j dp B dpidp 2 dp 3 x 

O-Pi . P2, pr, N q , N q )Rf{p B ,pu P2, P3;N q ,N q ), (23) 
and obtain 

N Mj (N qi ,N qt ) = 2 N qi q 2 K^ j>qi9 j M , (24) 

9 192 

N Bj (N qi ,N qi )= 2^**^.*^- (25) 

919293 

Summing over different species of mesons and those of the 
baryons, respectively, we obtain the average total numbers 
of mesons and baryons produced in the combination of the 
system of N q quarks and N q antiquarks as 

N M (N q , Nq, Va ) = N qq y M (N q ,N q , V* ), (26) 
N B (N q , N q ,^s) = N qqq y B {Nq, N q , ), (27) 

where N qq = N q N q and N qqq = N q (N q - l)(N q - 2) are 
the total number of qq pairs and that of qqq systems, re- 
spectively. The factors YM(N q ,N q , yfs) and y B (N q ,N q , yfs) 
represent the probability for a particular qq from the sys- 
tem consisting of N q quarks and antiquarks to combine 
with each other to form a meson and that for a qqq to form 
a baryon, respectively. We emphasize in particular that 
Eq. (26) does not mean that the average yield of mesons 
is proportional to the product of the number of quarks and 
that of antiquarks since the factor yM can depend strongly 
on N q and Af.j. This is because, for a given q, the larger 
N q and/or A^, the more possibilities for the q to combine 
with the others to form a hadron, and thus the smaller the 
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probability for it to combine with the given q to form the 
meson. In fact, it can in general be expected that y M should 
be more or less inversely proportional to N q and/or Nq and 
the final result for (Nm) should be roughly proportional to 
N q + Nq. A similar conclusion holds for {N B }- 

In terms of the total average numbers of mesons and 
baryons, the average number of a specified meson My and 
that of a specified baryon B y are given by 



flavor composition; and for baryons, 



N Mj (N qi ,N 9i ) = J] ^nfi (N qi ,N 9t ) 

xN M (N q ,Nq,^), 
N Bi (N qi ,Nq t )= Z^^3^.^) 

qiqiq, ^ 

xN B (Nq,N q ,^s). 



(28) 



(29) 



The parts describing the flavor dependence of the kernel 
functions, ^^(Nq.Ng,) mdK^^N^Nq,), have to 
guarantee flavor conservation in the combination process. 
Hence, they contain the Kronecker (Ts and constant factors 
Cm/ and C Bj . For example, for 7r + and p, they are given by 



Cr, 



Roiold + Ro/d) for J p = (l/2) + baryons, 
1/(1 + Ro/d) for J p = (3/2) + baryons, 



(34) 

except that C A = C 2 o = R /d/(1 + 2R /d), C 2 .o = 
1/(1 +2R 0/D ), and C A ++ = C A - = C n - = 1. Here, 
Ro/d stands for the ratio of the J p - (l/2) + octet to the 
J p = (3/2) + decuplet baryons of the same flavor compo- 
sition. The two parameters Ry/p and Ro/d can be deter- 
mined by using the data from different high-energy reac- 
tions [1, 13, 38]. 



C. Modeling P(N q . , (N q . > , yfs) 

We consider three flavors u, d, and s of quarks and anti- 
quarks. Inside the system of N q quarks and N q antiquarks, 
we suppose that each quark can take flavor u, d, or s with 
given probability p„, pd, or p s independent of the others. 
In this case, the numbers of u, d, and s quarks inside the 
system at a given N q obey the multinominal distribution, 
i.e., 



T ,<?192 
(/) 



— C pi& q\,i$ q2,i$ qi,d ^qi,u^q2,d^qs,u 
+ $qiA8q 2 ,uSq 3 ,u)- 



(30) 



We recall that, in the factorized case considered here, the 
flavor-dependent part 9?^ qiq2 {N qi ,Nq?) of the kernel func- 
tion represents the probability for the specified q\qi with 
the specified flavor q\ and qo from the system consisting 
of N q quarks and N q antiquarks to form the specified me- 
son Mj under the condition that they are known to form a 
meson. Although we can not prove it, it is very unlikely 
that this probability still depends strongly on the environ- 
ment. We therefore consider the simplified case where Cm s 
is taken as a constant independent of N q or N q . The same 
applies to qmqs (N qi ,N qi ). In this case, we have 

N M] (N q „N qi ) = C M ^-N M {Nq,Nq, Vj), (31) 



N, 



N B j(N qi ,N qi ) = N iter C Bj -jp^N B (N q ,N q , yfs), (32) 



where Nj ler stands for the number of possible iterations of 
9i?2^3 which is 1, 3, and 6 for three identical flavor, two 
different flavor, and three different flavor cases, respec- 
tively. 

In the case when only J p = 0~ and 1" mesons and J p = 

i + -i + 

| and | baryons are considered, we have, for mesons, 

_ f 1/(1 +Rv/p) for J p = (T mesons, 
M < ~ \ R V /p/(l + Rv/p) for J p = I - mesons, ( ' 

where Ry/p represents the ratio of the J p = 1~ vector 
mesons to the J p - 0" pseudoscalar mesons of the same 



B(Nq-N q ) = 



Nq\ 



N u \N d \N s \ 



N„ N d N, r 
Pu Pd Ps °N r N u +N d +N s , (JJj 



where p u = p d = 1/(2 + A q ), p s = A q /(2 + A q ), and A q 
is the effective strangeness suppression factor for quarks. 
Similarly, for the antiquarks, we have 



B(Nq-Nq) 



N, 



NnlNdlN, 



~Pa Pd Ps 6 N s ,Nr,+N 3 +Ni, (36) 



where p a — pj = 1/(2 + A), p- s = A/(2 + A), and A is the 
strangeness suppression factor for antiquark production. 

In general, in high-energy heavy-ion collisions, the sys- 
tem contains the contributions of the net quarks coming 
from the incident nuclei. Hence the effective strangeness 
suppression factor A q for the quarks is different from A for 
the antiquarks, which do not have influence from the net 
quarks. Here, we keep them as distinguished from each 
other so that we can apply the results to different cases. 
Furthermore, we emphasize that the system considered 
corresponds to a quark-antiquark system produced in an 
AA collision in a limited kinematic region. The system is 
supposed to be a small part of the whole quark-antiquark 
system produced in the reaction so that the influence from 
the global flavor compensation is considered to be neg- 
ligible. The global flavor compensation can have some 
influence on the flavor correlation in hadron production. 
Such a case was for example discussed in [5] for baryon- 
antibaryon flavor correlation in e + e annihilation where the 
number of quarks was of the order of tens and the aver- 
age yield of baryons in an event was less than one. It was 
found that, even in that case, the global flavor compensa- 
tion does have some effect on the flavor correlation but the 
effect is not very large. Hence, for simplicity and clarity, 
we neglect them in the discussion here. 
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At given N q and N q , we average over the distribution 
of the numbers of quarks and/or antiquarks for different 
flavors. It can easily be shown that, for q\ + qi + ^3, 

Y J N qi B(N qi ;N q ) = N qPqi , (37) 

J^N^N^BiN^Ng) = N qqqPqiPq2 p qy , (38) 

YjN qi (N q , - \)N qi B(N q -N q ) = N qqq p 2 qi p q2 , (39) 

and similarly for others, so we obtain 
N Mj (N q ,N q , V^)=C Mj p qi p q2 N M (N q ,N q , V?), (40) 
N Bj (N q ,N q , yR)=NiterC Bj p qiPq2 p q3 N B (N q ,N q , V^).(41) 

For a subsystem of quarks and antiquarks in a given 
kinematic region in AA collisions at given energy rfs, (N q ) 
and (N q ) are fixed while N q and N q follow the distribu- 
tions P q (N q ; (N q )) and P q (N q ; (N q )), respectively. Hence, 
we need to average over these distributions and obtain 

(N Mj )((N q ), (N q ), ^)=C Mj PqiPh (N M ), (42) 
(N Bj )((N q ), (N q ), ^)=N iter C Bj p qiPq2 p qi (N B ), (43) 

where (Nm) and (N B ) are functions of (N q ), (N q ) and -\fs 
and stand for the average total number of the mesons and 
that of the baryons produced in the combination process. 
They are given by 

(N M )=J]P q (N q ;(N q ))P q (N q -(N q ))N M (N q ,N q , V?), 
(N B )=J]P q (N q ;(N q ))P q (N q ;(N q ))N B (N q ,N q , 

We see that, in this case, for the directly produced 
hadrons, the ratios of the yields of different mesons, those 
of different baryons and those of the antibaryons separately 
are constants depending on the parameters A, A q , Ry/p and 
Ro/d- In general the effective strangeness suppression fac- 
tor A q for quarks contains the influence from the net quark 
contributions and can be dependent on (N q ) and (Nq). This 
leads to a dependence on (N q ) and (N q ) even for such kinds 
of hadron yield ratios. In the case that the net quark con- 
tribution is negligible, we have A q - A; these kinds of 
hadron yield ratios become constants independent of (N q ) 
and (N q ). This should be the case for a subsample in the 
central rapidity region of the bulk quark-antiquark system 
produced in AA collisions at very high energies such as 
those at the LHC. In this case, we have also (N B ) = (N B ) 
and this, together with p q . = p q ., leads to (N B .) = (Ng ). 
These are predictions that can be checked at the LHC. 

D. Including the decay contributions 

Including the decay contributions, we calculate the 
yields of different hadrons in the final state. We denote the 



decay contribution from a hadron /z, to hj by Br(hi — » hj) 
and obtain 

(N{) = (Nhj) + J] B < h i "» hjW hi ), (44) 

where we use the superscript / to denote the results for 
the final hadrons to differentiate them from those for the 
directly produced hadrons. Here we consider only the in- 
fluence from the decay of the short-lived hadrons but do 
not consider the influences from the final-state interactions 
of the hadrons. 

The value of Z?r(/j, — > hj) can be obtained easily from 
the materials given by the Particle Data Group [39]. In 
the following, we take the strong and the electromagnetic 
decays into account. For most of the hadrons, the results 
look very simple. In the case in which only J p = (T and 1~ 
mesons and J p = (l/2) + and (3/2) + baryons are included, 
the average yields of final hadrons, e.g., K + , p, and A, are 
given as 

(A^ + > = <Nk+) + |<Afro> + 1<Nk») + 0.489<AV>, (45) 
(A^> = (N„) + (N A „) + ^(7Y A+ ) + ^(N A o), (46) 

(N J A ) = (N A ) + (N^) 

+ 0.883<iV z .o> + 0.94«7Ve. + > + <A^-». (47) 

We consider the case discussed in Sec. II C and sub- 
stitute the results for (Nmj) and those for (N B ) given by 
Eqs. (42) and (43) into the above equations and obtain 

f 0A89R V /p 
(N J K+ ) = p u p,(Nm) + -: ^PsPs(N M ), (48) 

1 + Ky/p 

(N f p ) = 4pl(N B ), (49) 
f (530+12R O /d 5.64 \ , 

where we have taken p u = p^ and p r , — pj. 

Taking R V / P = 3 according to the spin counting, but 
Ro/d = 2 since decuplet baryon production is observed 
much more suppressed [1, 13, 38], we then have 

(N f Kt ) = (puPs + 037p s p s )(N M ), (51) 
(N{) = l.lAp 2 uPs (N B ). (52) 

In Table I, we show the results obtained for baryons 
and antibaryons and those for strange mesons in differ- 
ent cases. We see that, like for protons, the results for 
many final baryons look even simpler than those for the 
directly produced cases since the corresponding decuplet 
baryons decay strongly to these baryons. This makes the 
average final yields for these baryons independent of the 
ratio Ro/d- In Table I, we also present the results for the 
simple case without net quarks. The results obtained if we 
take Ry/p - 3 and Ro/d = 2 are also given. 
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TABLE I: Average yields of the directly produced hadrons and those including decay contributions. Here, in the second column, we 
show the results for the directly produced hadrons. In the third column, we show the results when the strong and the electromagnetic 
(S & EM) decay contributions are taken into account. The fourth column shows the results when the net quark influence is negligible. 



In the last two columns, we see the results for the case when R 



v/p 



3 and R, 



O/D 



Hadron 



Directly produced 



With S & EM decays 



N": 



(/L = X) 



^V/P 



3,R, 



O/D 



N net ^ 



N ne, 



3R QID _3/»r v 

2/A7 \ 

3R OID „3/AT.\ 
3R OID -3/VA 

i + Ro /d P«Wb) 

2/ AT \ 

ur^PuPsWb) 

3Rq/D 2/ AT \ 

TTr^PhPANb) 
p\{Ns) 



*Pl(N B ) 
4pl(N B ) 

3puP 2 AN B ) 
pKNb) 
4pl(N s ) 

^pI(n- b ) 

3puPl(N s ) 
2pup](N s ) 
p]{N B ) 



ft<M»> 



■{N B ) 



(2+-)) 

ft™ 
ft™ 



4^<JV B ) 
3p uP 2 s {N B ) 
3p u p 2 (N B ) 
pIWb) 

4pl(N B ) 

4pl(N B ) 
2puP 2 {N B ) 
2p«P 2 (N B ) 

pUNb) 



ww(*b) 

WwWb) 

ft<*'> 
ft<*»> 



(2+.1) 3 



m^PuPs<N M ) 

u^PuPs<N M ) 
u^PuPANm) 
uk^PuPANu) 

^- d pIps<n b ) 
^pIpANb) 
^- d pIps<n b ) 

3^o/D 2 /AT \ 

uS^PuPs(N B ) 

^epIMNb) 



p.a(i + ^M)<**> 

tts-P s P-ANm) 



(2+-)) ; 

(2+.t) : 

(2+,l)2 
,1 

(2+.t) : 



4 .0 + T^)<^> 



+ 



I+2J? 0/D l+Ro/D 
0.35 . 0.1 



, 1+2R /D 
(6Ro/D 



1+ ^)pIpAn b ) 
'-)p 2 ,pAn b ) 
;)pIPs(N b ) 
'-)pIpAn b ) 



l+Ra/D 



A /5.30+12K 
(2+,l)H 1+2R , 

i f 0.35 
(2+i)- 1 1 l+2R , D 

(2+/1) 3 ( 



1+Ko/fl 



1+Ko/fl 
6«o;n 0.36 



+ 



(2+A) 3 V 1+2K /D 1+Ko/fl 



l+«0/D 



(2+i) 3 ( 



0.35 



(2+A) 3 V l+2R ,D 



)p\pAN B ) 
)pIp-An-b) 
)pIp-ANb) 

I8+3* 0/d \ 2 /at \ 



, 5.64 

.18+3K 0/fl 
l+*0/D 
0.36 



l+«0/D 



5.64 
i+lo/D 



(2+.1) 3 V 1+2K /Z> l+Ro/D 
6Rq/d _,_ 0.36 



!+!l) 3 ( 



+ 



(2+i) 3 Vl+2Ro/D l+«o/D 



0.35 

(2+i) 3 V l+2^ 0/o 



(2+i) 3 ( 



+ 



)(N b ) 
'-)(Nb) 
)(N b ) 

)(N B ) 
'-)(N B ) 



p u p^+031X q ){N M ) 
p B p s (l + 031X)(N M ) 
p uPi {l+0.26A q )(N M ) 
p aPs {\+0.26X)(N M ) 
\p s p-AN M ) 
I.IApIpANb) 
2.U P 2 lPs (N B ) 
2,52pIpAN b ) 
2.UpIpANb) 
I.IApIpANb) 
2A3p 2 P AN B ) 
2,52pIpANb) 
2A3 P IpANb) 



(2+A) 



f(N M ) 



A+0.37A 2 / at \ 
A+OlbA 1 , », x 



-1+0.26.1- , N x 

!(afe)V*> 
££(Nb) 



(2+A)* 



j(N B ) 



We do not list the corresponding results for pions in Ta- 
ble I since the corresponding expressions are quite long. 
This is because the pion receives contributions from the 



decays of almost all the other mesons, baryons, and an- 
tibaryons. For example, for ji + , we have 



<A&> =<AW + <JV> + <A>> + |«iV*o> + (N K »)) + 0.2734<Ay + 0.9073<^) + 0.9274{N, r ) + 0.1568<iV,> 

+ 0.94«7Y 5: , + > + <iV £ , + » + 0.0585«^V z .o) + </V £ . » + ^«^V 3 .o) + <iV s , + » + <^V A++ ) + i«iV A+ > + </V A „» + (N A+ ). (53) 
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By inserting the results given by Eqs. (42) and (43), we obtain 

f 1.71+2.91/?™ 0.49 + 0.16/?™ 2 /?™ 

= 1 5 PuPa(N M ) + : PsPs(N M ) + 7 1 - — PuPs(N M ) 

1 + Ky/p 1 + Ky/p J I + Ky/p 

2 Ry/p 2.82 0.35 , 2.82 0.35 , 

+ ~ t—^—PuP- s (Nm) + (t—= + )pIp s (Nb) + ( — + -—^-)pIps(Nb) 

2 2/»t x 2 9 2 + Ro/d 3... . 2+Ro/d -x,. T . .... 

1 + "0/Z5 1 + «0/D 1 + "O/fl 1 + KoiD 

For systems without net quarks, we have p u = pa = pd = p d and p s = p^, and thus we obtain 

f 1.71+2.91/?™ 2 0.49 + 0.16/?™ 2 4 /?™ 

<A W = TTff + i"T~p p ~* {Nm) + ? i i p PuPANu) 

1 + Ky/p 1 + Ky/p J 1 + /vy/p 

5.64 0.70 , 4 , 4 + 2/?o/d , 

If we take R v /p = 3 and /?o/d = 2, we have 

<j\£> =2.6lp 2 u (N M ) + p B ^ s (iVM) + 0.24^</V M ) + ^p tt p 2 s (N B ) + 2.02/^ s </V fl > + ^p 3 u (N B ) 
2.61+A + 0.24A 2 8/3 + 2.02/1 + 4/3/1 2 

= (2 + , )2 + ( 2 + ,)3 W»>- ( 56 ) 



From these results, we see clearly that there exist 
many simple relationships between the yields of different 
hadrons. These are the characteristics for hadron produc- 
tion in the combination mechanism. We will list some of 
these simple relations in the following. Before doing that, 
we first discuss the net quark influences in the next section. 

E. Influence of the net quarks 

In a heavy-ion collision at high energy, the produced 
quark-antiquark system consisting of the newly produced 
quarks and antiquarks and the net quarks from the incident 
nuclei. For a subsample of this quark-antiquark system in 
a given kinematic region, we have, in general, 

(N g ) = (N 9 ) + (N'p. (57) 

Both the momentum and flavor distributions of these net 
quarks are different from those for the newly produced 
ones, and this leads to observable effects in the final 
hadrons produced in hadronization. We expect that they 
have influences on the following aspects: 

(i) The difference in momentum distribution leads to dif- 
ferent jM(N g ,Nq, -\fs) and jB(N q ,Nq, yfs), as seen clearly 
from Eqs. (22) and (23). Furthermore, since (N q ) > (Nq), 
the average number of baryons, (Nb), should be accord- 
ingly larger than (N B ). The ratio (N b )/(Nm) should in gen- 
eral depend on (N n q et ) I (N q ) . 

(ii) The distribution of the number of quarks, N q , at a 
given (N q ) is different from the corresponding distribution 
of the antiquarks. The distribution of the number of the 
net quarks, N q let , at a given (N' q wt ) is different from those 



for the newly produced quarks and/or antiquarks. This 
leads to a difference between the quark number distribu- 
tion Pq(N q , (Nq), -\fs) and the antiquark number distribu- 
tion P q (N-q, (Nq), js). 

(iii) The flavor distribution of quarks is different from 
that for the antiquarks. The net quarks take only two fla- 
vors, u and d. At a given N q et , the average numbers of u 
and d net quarks are determined by the numbers of protons 
and neutrons in the incident nuclei. For a given AA colli- 

— net — net 

sion, N u : N d = (A + Z) : (2A - Z), where A and Z are 
the numbers of nucleons and protons, respectively, in the 
incident nucleus A. The numbers N"" and N n A et follow a 
binominal distribution with p" u e ' : p'f = (A+Z) : (2A-Z). 
For the newly produced quarks or antiquarks, at a given 
pjnew or jy_^ tne nurm -,ers of them of different flavors fol- 
low the multinominal distribution as given by Eq. (35). 
Hence, including the net quark contribution, the distribu- 
tion of the numbers (N u , Nd, and A^) of the different fla- 
vors (u, d, and s) of quarks at a given number of quarks 
(Nq = N u + Nd + N s ) is also different from the correspond- 
ing distribution for the antiquarks. 

The detailed calculations of the influences of these ef- 
fects on the hadron yield ratios in the combination mech- 
anism depend on the particular models. In this paper, we 
present a rough estimate of, at least, the qualitative ten- 
dency of these effects by using the following two approxi- 
mations. 

First, we approximate that the flavor distribution of the 
number of quarks in the subsample of the system is inde- 
pendent of that for the antiquarks. That for the antiquarks 
is given by the multinominal given by Eq. (36). For the 
quarks, we approximate it by a multinominal distribution 
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in the same form as that for the antiquarks but with dif- 
ferent probabilities of the different flavors. For the newly 
produced quarks, the flavor distribution should be 1 : 1 : A 
for u, d, and s quarks, the same as those for the antiquarks, 
but for the net quarks, it should be (A + Z) : (2A - Z) : 0. 
We take the average and obtain, at given N q and N™', 



1 



Pu 



Pd 



2 + A 
1 

2 + A 
A 

2 + A 



(1- 



(1 



(1 



jwnet 

q 

N q 

N net 

1 > 

N, ' 



) + 



) + 



3A N q ' 

2 A-zK;" 

3A N a 



(58) 
(59) 
(60) 



In the case in which the u and d difference is not very large, 

— net — net 

we neglect it and consider the case in which N u : N d = 
1 : 1 . In this case, we have, at given N q and N'™', p u : pd ■ 

Ps 



= 1 : 1 : A„ and 



= ^=41 

N u 



N net 



TN a -Nl 



3 



(61) 



Under this approximation, we see that Eqs. (40) and (41), 
the results we obtained in Sec. II C for directly produced 
hadrons in the combination of a quark-antiquark system 
at given N q and Nq, are still valid. We only need to note 
that the p q . in this case is a function of N q and N n q et as 
given by Eq. (61). The influence of the isospin violation 
in net quarks can manifest itself in the difference between 
the average yields of hadrons belonging to the same charge 
multiplet. This can be studied separately in experiments. 
Since our purpose is a rough estimation of the net quark 
influence, we consider in the following first the simplified 
case where u and d are equal but leave the isospin differ- 
ence for future studies. 

Second, for a subsample in a given kinematic region 
of the bulk system produced in a given AA collision 
at given energy the averages (N q ) and {N q } are 
fixed and the numbers N q and N q follow the distribution 
P(N q ,N q ;(N q ),(N q ), -fi). The average numbers of the 
hadrons produced should be the average over this distribu- 
tion. In general such averages depend on the precise form 
of P(N q ,N q ; (N q ), (N q ), -*fs). In the rough estimations we 
made here, we approximate these averages by taking the 
corresponding values of the quantities at the averages (N q ) 
and (Nq), i.e., 

(N hj (N q , N q , VS)> * N hj ((N q ), (Nq), V?). (62) 

Under these two approximations, all the results pre- 
sented in the last four sections where we distinguish be- 
tween pq. and p q . apply and we can use them to make esti- 
mates of the effects of the net quarks. 



F. Ratios of the yields of different hadrons 

From the results given in Table I, we see that there 
are many simple relations between the yields of different 



hadrons. In particular, for the directly produced hadrons, 
such relationships are very simple. Even for the final 
hadrons, although the decay influences are often very 
large, there still exists a set of simple relations between 
them. For example, independent of the values of Ry/p 
and Ro/d, for the final hadrons where contributions from 
strong and electromagnetic decays are taken into account, 
we have 



(N f p ) 
(N f p ) 

<^ 
<> 

(N{) 



(n{,) 



= -At 



4* 



(Ho) 



<N f p ) 



= -A 



-A 



_ , 2 + A q ,3(N B ) 
A ' 2 + A > (NbY 



(%) 

<^A> 

<<> 

(N f n .) 



2 + A' Aq (Nb) ' 

_ , 2 + A q ,3,Aji(N B ) 
'2 + A> ^A q > (N B )' 

(2 + A q \3{ A ^ (N B ) 



^2 + A'^Aa' 



(63) 
(64) 
(65) 
(66) 
(67) 
(68) 
(69) 
(70) 



In the case in which net quark contribution is negligible, 
we have, A q — A and (N B ) = (N B ), so that 



(nL) 

(N f p ) 

<-> 
(N f p ) 



(n1 ) 

(N f p ) ' 



-A 1 



-A' 



(71) 



(72) 



and (N f s )l(N f B ) = 1 for all the different types of Bj. 

In the case in which Ry/p - 3 and Ro/d = 2, we have 
more such simple relations such as 



(N^ 
(N f p ) 

<<> 



/j,w 1 +0.37,1 x 
' A " 1 +0.37iJ' 



= 1.935^1, 



= 1.935X, 



(73) 



(74) 



(75) 
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These relations are intrinsic properties of the combination 
models in the sense that they do not depend on the de- 
tails of particular combination models but are determined 
mainly by the basic ideas of the combination mechanism. 
They can be used to determine the free parameters and/or 
to test the mechanism. We note in particular the following 
two features. 

(i) For a system with equal average numbers of quarks 
and antiquarks, i.e., where net quark contributions are neg- 
ligible, these ratios are quite simple and can be tested by 
the data in extremely high energy collisions, e.g., at the 
LHC. 

(ii) With net quark contributions, most of these yield ra- 
tios of different hadrons depend on (N^ et )/(N q ). We see in 
particular that the ratios of the yields of hadrons to those of 
the corresponding antihadrons deviate from one in general. 
They tend to one for reactions at very high energies where 
the net quark contribution tends to vanish. This leads to an 
energy dependence of such ratios, even for particles such 
as fi and £2 + . Such a property for the combination mech- 
anism is different from what one expects from fragmenta- 
tion and can be used as a check to differentiate the different 
hadronization mechanisms. 

We can also build some combinations of the average 
yields of the hadrons and obtain simple results for some 
more sophisticated ratios such as, 



A = 



B = 



(NkWk-XN,,) 

<N„)(N K+ XNp) 
<N A )(N K -)(Nz + ) 



(Nfi)(N K+ )(N s -) 



= 1 , 



= 1 , 



(76) 



(77) 



(NpKNa-) 



(79) 



We see that they all lead to simple results and these re- 
lations are not influenced by the resonance decays except 
<p — > K + K , which slightly changes A and B from unity. 
For the case in which N' q K1 = 0, all four ratios are equal to 
unity. 

As a brief summary, we emphasize once more that the 
model that we consider in this section is intended to be a 
general case based on the basic ideas of the combination 
mechanism. The purpose is to concentrate on the hadron 
yield correlations in the combination models. No effort 
is made to study other properties such as momentum dis- 
tribution, etc. The results obtained follow from the ba- 
sic ideas of the combination mechanism and a number of 
assumptions, simplifications, and/or approximations such 
as the factorization of flavor and momentum dependence 
of the kernel functions, the flavor independence of the 
quark-antiquark momentum distribution, the independent 
production of different flavors of quarks and antiquarks, 
and the approximations made in considering the net quark 
influences. These results do not depend on the detailed 
form of the momentum dependence of kernel functions 



and/or the momentum distributions of the quarks and an- 
tiquarks. They even do not depend on whether the quark 
number conservation or depletion is imposed in the combi- 
nation process. Such a conservation or depletion of quark 
number influences the relationship between the average 
number of mesons (or baryons, or antibaryons) produced 
and the number of quarks and/or antiquarks participating 
in the combination process (see in particular the discus- 
sion in [40]) but does not influence the hadron yield ra- 
tios if the factorization is assumed. These results should 
also be valid in the combination models discussed in the 
literature wherever these assumptions and/or approxima- 
tions are also made (explicitly or implicitly). Some of 
them should even be common in these different models 
[9-13, 23]. In fact, some of the relationships presented 
above have also been derived in these literature [9-13, 23]. 
For example, relations similar to those given in Eqs. (76)- 
(79) have been obtained in [9]. These relations can be used 
to test the validity of the combination mechanism and the 
assumptions made. We will compare them with the data 
available in the next section. 



III. COMPARISON WITH DATA 

There are already quite abundant data available from ex- 
periments in a quite broad energy region, from low SPS 
energies to RHIC and LHC energies [18, 21, 22, 24-36]. 
We compare the results obtained in the last section with 
these data in the following. 

As mentioned earlier, the study presented in Sec. II is 
intended to be a general case for the combination of a 
quark-antiquark system consisting of N q quarks and Nq 
antiquarks. No effort is made to ascertain whether the com- 
bination mechanism dominates the production of hadrons 
in the given kinematic region in AA collisions. In this sec- 
tion, we choose the data in the central rapidity regions in 
different AA collisions at different energies. The agree- 
ment and/or disagreement of the theoretical results with 
these data should give us a signature of whether the combi- 
nation mechanism with the above-mentioned assumptions 
and/or approximations is applicable. Also, since some of 
the theoretical results depend on more inputs and some of 
them depend on fewer inputs, we make the comparison at 
different levels. 



A. Comparison with LHC data 

At the first level, we consider a subsample of the quark- 
antiquark system in the central rapidity region produced in 
AA collisions at very high energies. We suppose the ener- 
gies are very high and the subsample that we consider is 
only a small part of the whole quark-antiquark system pro- 
duced in the collision process so that the influence of the 
net quarks and that from the global flavor compensation 
are negligible. In this limiting case, the results for the ra- 
tios of the yields of different hadrons are divided into three 
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classes. 

In the first class, we consider the ratios of the yields of 
hadrons to those of the corresponding antihadrons. Such 
ratios are unity, independent of any parameter. This can 
be considered as a criterion for the validity of this limiting 
case. Results from LHC experiments can be considered as 
an example for this case. In the first three lines of Table II, 
we show the available experimental results for the parti- 
cle to antiparticle ratios such as (A^_)/(A^+), (A^-)/(A^ + >, 
and (Np)/(Np) at mid-rapidity. The data are obtained from 
Ref. [21]. We see that they are indeed very close to unity. 

TABLE II: Hadron yield ratios obtained for N" a = compared 
with data from the LHC in Pb + Pb collisions ai ^ = 2.76 TeV. 
The data are taken from Refs. [21, 22]. The experimental result 
for K~ (ic is used to determine the strangeness suppression factor 
A. 



Ratios 


Data 




Calculations 


7T~ /tt + 


1.000 + 


0.080 


1 


K~/K* 


0.987 ± 


0.076 


1 


PIP 


0.995 + 


0.077 


1 


<f>/K + 






0.278 


K°/K + 






0.959 


A/p 






0.832 


E-/p 






0.139 


Or/p 






0.020 


K-jn 


0.155 ± 


0.012 


0.155 


p/n + 


0.045 ± 


0.004 


0.043 


A/7T+ 






0.036 


E-/7T+ 


0.005 ± 


0.001 


0.006 


Q.-/n + 


0.001 J 


: 0.0002 


0.001 


p/K + 






0.275 


A/K + 






0.229 


Z-/K+ 






0.038 








0.005 



In the second class, we consider the ratios of the average 
yields of hadrons such as (p/K + , A/p, and so on as shown 
in the second part of Table II (from the forth line to the 
eighth line). Because in this case we have (Nb) = (N B ) 
and A g = A, these ratios depend only on one free parame- 
ter, the strangeness suppression factor A, and are indepen- 
dent of the particular models. We can fix A by using the 
data for one ratio and make predictions for other particle 
ratios. Such results can be used to check the validity of the 
combination picture. 

In the third class, we consider the ratio of a specified me- 
son to a specified baryon. To obtain the results for such ra- 
tios, we need the input for (Nb)/{Nm)- This can be slightly 
different in different combination models. As an example, 
in the third part of Table II (from the ninth line to the end), 
we show the results obtained by taking 

(N B )/(N M ) = 1/12. (80) 
This is obtained by parametrizing the results for large 



N q - N q of the Monte Carlo generator (SDQCM) based 
on the combination rule developed by the Shandong group 
[13, 15, 37] which has reproduced the data well. This 
parametrization is valid with high accuracy for N q larger 
than, say, 100. The parameter A is taken as A — 0.43 by fit- 
ting the data of K~/n~. We see that the theoretical results 
agree with the LHC data whenever available. 



B. Comparison with RHIC and SPS data 

At the second level of comparison, we consider the case 
where N q et + 0. In this case, we need a further input of 
(N q e ')/(N q ) which describes the strength of the net quark 
influence. Clearly, this ratio depends on the type and the 
energy of the incident nuclei, and also on the kinematic re- 
gion that we consider. For example, to make a good com- 
parison with the data from the RHIC and the SPS [18, 24- 
36], we need to take such effects into account. 

In practice, to carry out the calculations, we first fix A q 
and A by using the data for (Ns-)/(N p ) to determine A q and 
(N§+)/{Np) to determine A, and we then derive (N q e ')/(N q ) 
using Eq. (61). In Table III, we show the results obtained 
by taking the data for dN/dy at y = (where y denotes 
the rapidity of the hadron). These results can be used to 
calculate the hadron ratios in the mid-rapidity regions. 



TABLE III: The fixed A and A q and the derived (N™')/{N q ) at 
RHIC and SPS energies. 



Energy 


RHIC (GeV) 


SPS E} )ewn 


(A GeV) 


200 


130 


62.4 


158 


80 


40 


30 


A 


0.425 


0.412 


0.369 


0.499 


0.567 


0.540 


0.645 


** 


0.397 


0.375 


0.328 


0.255 


0.232 


0.193 


0.192 




0.056 


0.074 


0.095 


0.434 


0.529 


0.587 


0.640 



From the table, we see that {N q e ')/{N q } is small at RHIC 
energies but quite large at SPS energies and the expected 
effects should be large at those energies. 

Using the obtained A and A q values as inputs, we cal- 
culate the hadron yield ratios that are independent of 
(Nb)/(Nm)- The results are given in Table IV. The cor- 
responding data are from Refs. [18, 24-35]. 

We see in particular that the ratio of the yield of hadrons 

to that of the corresponding antihadrons is not unity in this 

f f 

case. We also note that the ratio such as (N K -)/{N K+ ) = 
{A q /A)(l + 0.37/l)/(l + 031 Ag) is a good example to show 
the change of the effective strange suppression for quarks. 
The results should decrease monotonically with increasing 
{N n q e ')l{N q ). There are data available for {N f K .) I {N f K +) at 
different energies [18, 24, 25] and the data show clearly 
that the ratio increases with increasing energy. This qual- 
itative tendency is consistent with the effect of net quark 
contribution since the relative influence of the net quarks 
becomes smaller at higher energies. 

To calculate other ratios, we need (N b )/(Nm) and 
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TABLE IV: Comparison of the calculated results for hadron yield ratios that are independent of (N b )/{Nm) with data from RHIC and 
SPS energies [18, 24-35]. The experimental result for E~/p and that for E + //> are used to determine the strangeness suppression factor 
A q and A, respectively. 





V? (GeV) for Au + Au at RHIC 


Ebeam (A GeV) for Pb + Pb at SPS 


Ratio 


200 


130 


62.4 


158 


80 


40 


30 




data 


theory 


data 


theory 


data 


theory 


data 


theory 


data 


theory 


data 


theory 


data 


theory 


Rip 


0.94 ±0.15 


0.82 


0.93 ± 0.34 


0.80 


0.82 ±0.18 


0.71 


0.98 ± 0.23 


0.97 


1.22 ±0.30 


1.10 


1.31 ±0.35 


1.05 


1.31 ±0.41 


1.25 




0.14 ±0.03 


0.14 


0.13 ±0.03 


0.13 


0.10 ±0.03 


0.10 


0.19 ±0.05 


0.19 


0.24 ± 0.07 


0.24 


0.22 ± 0.07 


0.22 


0.31 ±0.15 


0.31 


Ci + /p 




0.019 




0.017 


0.017 ± 0.005 


0.013 


0.042 ± 0.020 


0.031 




0.046 




0.039 




0.067 


Np 


0.91 ±0.15 


0.77 


0.90 ± 0.30 


0.73 


0.78 ±0.16 


0.63 


0.37 ± 0.09 


0.49 


0.45 ± 0.08 


0.45 


0.37 ± 0.06 


0.37 


0.35 ± 0.06 


0.37 


sr/P 


0.12 ±0.02 


0.12 


0.11 ±0.03 


0.11 


0.08 ± 0.02 


0.08 


0.05 ± 0.01 


0.05 


0.04 ± 0.01 


0.04 


0.03 ± 0.01 


0.03 


0.03 ± 0.01 


0.03 


Qr/p 




0.016 




0.013 


0.010 ± 0.003 


0.009 


0.005 ± 0.001 


0.004 




0.003 




0.002 




0.002 


K-/K+ 


0.96 ±0.17 


0.94 


0.92 ± 0.09 


0.92 


0.86 ± 0.09 


0.90 


0.57 ± 0.05 


0.55 


0.48 ± 0.04 


0.46 


0.38 ± 0.04 


0.40 


0.37 ± 0.04 


0.34 


K°/K° 




0.94 




0.92 




0.90 




0.54 




0.44 




0.39 




0.33 


K°/K + 




0.96 




0.96 




0.97 




0.97 




0.98 




0.98 




0.98 



models. We also see that the antibaryon to baryon ratios in- 
crease with increasing energy, indicating that the net quark 
influence becomes smaller at higher energies. 

It was considered as a surprise that the net quark also 

f f 

influences <A^ + )/(A^_> significantly, as shown in Fig. 1, 
although the Q hyperon does not consist of u or d quarks. 
We see that from the lowest SPS to the highest RHIC en- 
ergy, {N^)/(Nq_) increases from about 0.4 to unity. This 
cannot be understood in the fragmentation models but can 
be naturally explained in the framework of quark combi- 
nation models. Because in the combination models more 
net quarks imply more chances for the antistrange quarks 
to meet quarks to form mesons, there will be more antis- 
trange quarks exhausted to form kaons than strange quarks. 
The probability for a strange antiquark to combine with 
other antiquarks to form an antibaryon is smaller than that 
for a strange quark to combine with other quarks to form a 
baryon. This leads to fewer O's than £Ts for the same num- 
ber of strange quarks and antistrange quarks. This quali- 
tative tendency is consistent with data, and from the figure 
we see also the quantitative results agree well with the data. 

We also compare the experimental results with the pre- 
dictions shown in Eqs. (76)-(79) for the "more sophisti- 
cated" ratios A = B = 1 and d~ = d 1 ^. From the data avail- 
able [24, 26-28, 36], we calculate these ratios and show 
the results in Table VI and Fig. 2 for d^ and d P a and A and 
B, respectively. We see that the data are consistent with 
d~ = d P a . The results for A and B at RHIC energies are 
consistent with unity while the error bars for those at SPS 
energies are too large to make a judgment. 



(Nb)I(Nm)- For this purpose, we parametrize {Nb)/{Nm) 



for the case where N" el * using SDQCM [13, 15, 37], 



and we obtain 

We found that this parametrization is a good approxima- 
tion to the results obtained from SDQCM for different N q 
in the range of N q > 100. In this model, quark number 
conservation or depletion is imposed in the combination 
process so that 3(N B ) + (N M ) = (N c/ ) and 3<A^> + (N M ) = 
(N q ) - (N'^'}, and we derive 

(N M ) { (N m ) 3 (N q ) ,/{ (N q ) I 

With Eqs. (81) and (82), we calculate the ratios of the 
yields of different baryons to mesons. The results are 
shown in Table V. The data are calculated from dN/dy at 
y — for different energies from Refs. [18, 24-35]. 

With the values of A and A q given in Table III and 
{Nb)/{Nb) derived from Eqs. (81) and (82), we calculate in 
particular the ratios of the antibaryons to the correspond- 
ing baryons at different energies. The results obtained 
are shown in Fig. 1 with the open symbols connected 
by different lines to guide the eye. The filled symbols 
with error bars are the experimental data taken from Refs. 
[24, 26-28, 36]. From Fig. 1, we see in particular that 
(N( + )/{N^.) > (N f 2+ )/{N^) > {N{)/{N f A ) > (N f p )/{N f p ). 
This was first pointed out by the NA49 Collaboration and 
was regarded as a distinct hierarchy of the antibaryon to 
baryon ratios [27]. The hierarchy can be naturally repro- 
duced by Eqs. (67)-(70) in the simple quark combination 
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TABLE V: Comparison of the calculated results for the baryon-to-meson ratios that are dependent of (Nb)/{Nm) with data from RHIC 
and SPS energies [18, 24-35]. 







Vs (GeV) for Au 


+ Au at RHIC 




Ebemn (A GeV) for Pb + Pb at SPS 


Ratio 


200 


130 


62.4 


158 


80 


40 


30 




data 


theory 


data 


theory 


data 


theory 


data 


theory 


data 


theory 


data 


theory 


data 


theory 


p/K~ 


O.Z / ± O.Oj 


0.25 


0.32 ± 0.0 / 


0.25 


0.31 ± O.Oj 


0.28 


U. 10 ± U.U2 


0.08 


0.0 / ± 0.01 


0.05 


0.04 ± 0.01 


0.04 


U.U2 ± U.U1 


0.03 


A/K~ 


0.26 ± 0.04 


0.21 


0.30 ± 0.09 


0.20 


0.26 ± 0.04 


0.20 


0.10 ± 0.01 


0.08 


0.09 ± 0.02 


0.05 


0.06 ± 0.01 


0.04 


0.03 ± 0.01 


0.03 


2+/K- 


0.04 ± 0.01 


0.03 


0.04 ± 0.01 


0.03 


0.03 ± 0.01 


0.03 


0.018 ±0.004 


0.015 


0.018 ±0.004 


0.012 


0.009 ± 0.003 


0.009 


0.006 ± 0.003 


0.008 


Cl + /K- 




0.005 




0.004 


0.005 ± 0.001 


0.004 


0.004 ± 0.002 


0.003 




0.002 




0.002 




0.002 


P/K + 


0.36 ± 0.07 


0.33 


0.42 ± 0.09 


0.36 


0.54 ± 0.08 


0.44 


1.00 ±0.15 


1.03 


1.22 ±0.18 


1.32 


2.05 ± 0.29 


1.79 


1.99 ±0.36 


1.93 


A/K + 


0.33 ± 0.05 


0.25 


0.37 ±0.10 


0.26 


0.42 ± 0.08 


0.28 


0.37 ± 0.08 


0.51 


0.55 ± 0.08 


0.60 


0.76 ± 0.09 


0.67 


0.69 ±0.10 


0.72 


Z-/K+ 


0.04 ± 0.01 


0.04 


0.04 ± 0.01 


0.04 


0.04 ± 0.01 


0.04 


0.049 ± 0.009 


0.050 


0.050 ±0.011 


0.054 


0.057 ±0.013 


0.050 


0.055 ± 0.014 


0.054 


Or IK* 




0.005 




0.005 


0.006 ± 0.001 


0.004 


0.005 ± 0.001 


0.004 




0.004 




0.003 




0.003 




VT NN (GeV) 



FIG. 1 : (Color online)The yield ratios of antibaryons to baryons 
at mid-rapidity at different energies. The filled symbols with error 
bars are the experimental results taken from Refs. [24, 26-28, 
36]. The corresponding open symbols represent the calculated 
results where the lines that connect these symbols are just used 
to guide the eyes. 



IV. SUMMARY 



TABLE VI: The deduced values of 4 and d" n at RHIC and SPS 
energies. Errors shown are total errors. The data are from Refs 
[24, 26-28, 36]. 



Energy 



200 GeV 
130 GeV 
62.4 GeV 
158A GeV 
40A GeV 




0.64 ± 0.07 
0.59 ± 0.08 
0.34 ± 0.04 
0.027 ± 0.009 



0.73 ± 0.23 
0.67 + 0.15 
0.32 ±0.14 
0.029 ± 0.016 



0.0013 ± 0.0005 0.0032 ± 0.0023 



(a) 



10 10 

VT NN (GeV) 



1.5 



00 1 



0.5 



(b) 



10 10 

VT NN (GeV) 



FIG. 2: (Color online)The correlation quantities A (a) and B (b) 
as the function of the collision energy. The experimental data, 
filled circles with total statistical and systematic errors, are from 
Refs. [18, 24-28, 36]. 



We study the hadron yield correlations in the combi- 
nation models. With the basic ideas of the combination 
mechanism and a few simplifications and/or assumptions 
based on symmetry and general principles, we show that 
the hadron yield ratios can be calculated and have a series 
of regular properties. These ratios are properties of the 



combination mechanism under these assumptions and/or 
approximations such as the factorization of flavor and mo- 
mentum dependence of the kernel function, the flavor inde- 
pendence of the momentum distributions, and the approx- 
imations made for the net quark contributions. They are 
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independent of the particular models where usually partic- 
ular assumptions are made for the kernel functions. A sys- 
tematic study of these ratios should provide good hints as 
to whether the combination mechanism is at work. Com- 
parisons with available data are made and predictions for 
future experiments are given. 
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